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Abstract 

V 

In  order  to  study  seismic  wave  propagation  through  laterally  varying 
sea  floor  structures,  a  software  package  has  been  created  to  generate  syn¬ 
thetic  seismograms  by  finite  differences.  The  elastic  wave  equation  can  be 
solved  in  two  dimensions  either  for  point  sources  in  cylindrical  coordinates 
or  for  line  sources  in  rectangular  coordinates.  Vertical  and  radial  varia¬ 
tions  of  the  elastic  parameters  are  allowed.' 

>The  package  includes  four  programs.  Input  to  the  system  consists  of  a 
short  file  containing  parameter  values  to  describe  the  model.  The  first  pro¬ 
gram  is  used  to  initialize  the  system  for  the  particular  model  being  used. 

The  source  arrays  and  velocity  matrices  are  each  computed  by  a  separate  pro¬ 
gram.  The  final  program,  which  actually  carries  out  the  finite  difference 
calculations,  includes  six  subroutines  to  implement  different  options  based  on 
alternative  finite  difference  formulations.  Two  different  kinds  of  output 
files  are  created  by  this  program:  one  or  more  snap-shot  files,  and  one  time 
series  file,  which  will  usually  include  more  than  one  series.  _ 
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Introduction 


This  report  describes  a  software  package  used  to  generate  synthetic 
seismograms  by  calculating  the  elastic  wave  equation  in  two  dimensions  by 
finite  differences.  The  original  finite  difference  code,  which  includes 
Kelly,  Alterman,  II an  and  Stephen  formulations  for  cylindrical  and  rectan¬ 
gular  co-ordinates  with  two  different  source  configurations,  was  developed 
and  written  by  R.  A.  S.  between  1978  and  1983.  The  absorbing  boundary 
subroutine  was  written  by  Mr.  S.  Emerman  and  is  reported  in  Emerman  and 
Stephen,  1983.  Applications  of  the  code  and  comparison  of  the  results  with 
the  reflectivity  method  are  discussed  in  Stephen  (1983).  L.  6.  assisted 
throughout  the  development  by  providing  programming  advice,  writing 
input/output  code  and  file  handling  procedures.  L.  G.  also  ran  or  directed 
a  series  of  'bench  mark*  tests  for  performance  evaluation  on  the  VAX,  on  the 
VAX  with  an  array  processor,  and  on  the  Cyber  205.  The  results  of  these 
tests  are  reported  in  the  Appendix.  M.  M.  H,  commencing  in  the  spring  of 
1983,  has  further  modified  the  code  to  be  more  'user  friendly'  and  to  be 
more  easily  transportable.  M.  M.  H.  also  prepared  this  manual  which  is 
being  written  in  order  to  facilitate  the  use  of  the  software  package  by  new 
users. 


The  program  development  is  an  on-going  process;  further  enhancements 
are  continually  being  made.  The  version  documented  and  released  here  is  a 
temporary  stoppi ng-pl ace. 

The  program  uses  a  variety  of  explicit  finite  difference  techniques  to 
solve  the  following  equation: 

pu*  =  (x+u)V(V-uWv2U  +  Vx(V-u)+Vu  X  (VXu)+  2  (Vp*V)u 

where  u  is  the  particle  displacement  vector 

x,u  are  Lame's  parameters,  and 
p  is  density. 


The  equation  is  solved  in  two  dimensions  either  for  point  sources  in 
cylindrical  co-ordinates  or  for  line  sources  in  rectangular  co-ordinates. 
Vertical  and  radial  variations  of  the  elastic  parameters  are  allowed.  The 
methods  have  been  compared  for  accuracy  with  the  reflectivity  method  by 
Stephen  (1983). 


II.  The  Model 


The  software  depends  heavily  on  the  definition  of  the  model,  which  is 
specified  for  the  computer  in  the  form  of  a  model  parameter  file.  A  diagram 
of  the  model  is  given  in  Figure  1. 

The  model  consists  of  three  horizontal  layers;  the  top  layer  is  water, 
the  middle  layer  is  the  sea-floor  boundary  layer,  and  the  bottom  layer  is  the 
sub-surface  stratum.  The  top  and  bottom  layers  are  assumed  to  be  homogeneous 
with  a  constant  velocity.  The  boundary  layer  has  a  varying  structure;  the 
velocities  within  this  layer  are  calculated  by  Program  FDBNY.  A  grid  is 
superimposed  on  the  entire  area  of  the  model.  The  grid  is  divided  into  NN 
depth  divisions,  and  MM  distance  (or  range)  divisions.  The  energy  source 
is  assumed  to  be  located  on  the  left  edge  of  the  model.  The  right  edge  and 
the  bottom  edge  of  the  model  are  either  rigid  or  absorbing  boundaries. 

There  are  two  categories  of  parameters  in  the  model  parameter  input 
file.  Some  parameters  describe  the  physical  model,  and  others  describe  the 
desired  processing  and  output.  A  description  of  all  parameters,  in 
alphabetical  order,  is  given  in  Table  I.  All  times  are  in  seconds,  and 
distances  are  in  kilometers. 
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Figure  1 

DIAGRAM  OF  MODEL 


RANGE-DEFTH  CROSS  SECTION:  IS0RB=0,  IFLAT=0 
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Table  I  Description  of  variables  in  model  parameter  file 
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III.  The  Programs 

The  package  includes  four  programs.  Input  to  the  system  consists  of  a 
short  file  containing  parameter  values  to  describe  the  model,  and  to  specify 
the  amount  of  output  desired.  The  final  output  includes  two  different  kinds 
of  files:  one  or  more  snap-shot  files,  and  one  time  series  file,  which  will 
usually  include  more  than  one  series. 

In  an  effort  to  minimize  the  computer  resources  used  for  a  given  model, 
minimum  array  sizes  are  computed  and  the  specifications  are  stored  in  COMMON 
files.  These  files  are  then  accessed  during  compilation  by  using  INCLUDE 
statements  in  the  programs.  The  charges  incurred  by  running  a  model  with 
larger  than  necessary  arrays  are  usually  greater  than  COMMLE  and  LINK  charges. 

The  first  program  in  the  sequence,  FDPREP,  checks  the  model  parameters 
for  validity,  and  determines  the  correct  dimensions  for  the  arrays  needed  by 
the  other  programs.  Two  files  are  created  with  COMMON  specifications 
correctly  formatted  for  INCLUDE  statements.  Descriptions  of  the  arrays  in 
these  COMMON  blocks  are  given  in  Table  II. 

The  second  program,  FDSORS,  calculates  sources.  The  program  creates  a 
file  containing  the  model  parameters  which  were  used,  and  the  source  arrays 
computed. 

The  third  program,  a  version  of  FOUNT,  finds  the  matrices  of  P-  and  S- 
velocity  squared,  and  the  density  matrix  of  the  boundary  layer  for  the  model. 
This  program  also  creates  a  file,  containing  model  parameters  used  and  the 
three  matrices  which  were  calculated. 

The  final  program,  FINDIF,  uses  the  model  parameter  file,  and  the  files 
produced  by  the  other  two  programs,  to  create  the  final  output.  The  snap-shot 
files  give  values  of  vertical  displacement  over  the  entire  grid  at  the 
specified  times.  The  time  series  give  values  of  vertical  displacement, 
horizontal  displacement,  or  pressure  at  each  of  the  selected  grid  points  for 
the  specified  time  interval. 

Program  FDSORS  must  be  compiled  using  the  specifications  for  COMMON 
block  COMSOR  which  were  found  by  FDPREP.  Programs  FOBNY  and  FINDIF  must  be 
compiled  using  the  specifications  for  COMMON  block  C0MFD8.  All  programs 
create  a  log  file  to  give  a  permanent  record  of  what  was  done.  A  diagram  of 
the  system  is  given  in  Figure  2. 
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Table  II 

Contents  of  COMMON  block  /C0MFD8/ 


DIMENSIONS  DEFINITION 

(MH»NN,2> .  Horizontal  <11  sol  acement 

(MNiNN,J)  Vertical  dlsolacement 

(MSW,NSW,2)  Used  (n  computation, 

<MSW,NSW,2»  only  w non  IFLAT  •  0 


(MM, 2)  or  (hhi 
( mm, 2  I  or  (mmj 

(M«| . . . 

(MM1 

•  MM , 2  I 

(MM) 


- PJGHTA 

RIGHTO 

- P/GA-I 

«MG3*l 

- POT  A 

nnre 

- of>T  »**1 

PnTB'*l 

AOATA 

SOC  S  A 

- 5nv  so-- 

- VP32 

V532 

- Bn? - 


These  arrays  are  used  In  the  comouta 1 1  on, 
only  when  ITRAN  «  1. 

When  IRECT  *  1,  the  dimensions  for  C  and  0 
are  (MM, 21!  otherwise  (MM  I . 

Used  In  computation, 

only  *n»n  IT&AN  >*  1  and  V51  *  0. 

Used  In  computation, 
only  when  IS0RI3  *  1 


( NMT , N  NT  «NKT I 
( MD I M ) 

( MO  I M )  .  - 


( MM,NPNDY» 
(MM,N3N0YI 
-(MM,NRNOY  I  - 


Time  series  array;  NMT* ( M0UT6N— MOUTST ) /M I  NT* 1 ,  etc 

Source  arrays  -  mdim  «  mm  |f  IFLAT  »  1? 

---  else  MDIM  *  4*m$w  ♦  2*NSW  -  8 

P-wave  velocity  snuared  In  boundary  layer 
S-waye  velocity  squared  In  boundary  layer 
Censlty  In  boundary  layer 

MINDY  -  NB— NA ♦ 2  IT  ITRAN  -  1,2; 

- else  NBNDY  *  NR-NAO.  -  -  ---  - 


Contents  of  COMMON  block  /C0I5S0R/ 


- RFR  IES - l2**NS0Rr.E  ) 

- rum* - ikoi",  mdim)  — 

nu"2  (kdim.moim) 


Comple*  array 

- Source  arrays  —  KOIM  *  K K  —  KSTRT  ♦  l 

MCI M  -  MM  If  IFLAT  •  1 

_  else  MO  I M  ■  4*MSW  ♦  2*NSW  -  ft 
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IV .  The  Files 


A.  Model  Parameter  File 


This  is  a  sequential  formatted  file,  most  of  which 
list-directed  READ  statements  (free-field  input). 
number  is  55.  The  record  contents  follow. 


is  read  with 
The  Fortran  unit 


1.  Record  1  consists  of  up  to  80  characters  containing 
identification  for  the  model. 


2.  Record  2  contains,  in  the  first  5  characters,  identification 
used  in  the  names  of  input  and  output  files. 

3.  Record  3  (free-field)  contains  values  of: 

IRECT,  IDENS,  I  FLAT,  IKELLY,  IEXPL,  IVERT,  ITRAN,  ISORB 

4.  Record  4  (free-field)  contains  values  of: 


MM,  NN,  KK,  KSTRT,  DELT,  DELR,  OELZ 

5.  Record  5  (free-field)  contains  values  of: 

KOUTST,  KOUTEN,  MOUTST,  MOUTEN,  NOUTST,  NOUTEN 

6.  Record  6  (free-field)  contains  values  of: 

KINC,  MINC,  NINC,  KMARK,  MI NC 

7.  Record  7  (free-field)  contains  values  of: 

VP1,  V SI ,  ROI ,  VP2,  VS2 ,  R02,  NA,  NB,  VPT,  VST,  ROT 

8.  Record  8  (free-field)  contains  values  of: 

NO,  MSW,  NSW 

9.  Record  9  (free-field)  contains  values  of: 

NSORCE ,  PLSWID,  TSWAVE 


B.  Boundary  velocity  fi 1 e 


This  file  is  created  by  FDBNY,  and  is  read  by  FINDIF.  It  is  a 
sequential,  unformatted  file,  created  with  binary  WRITE 
statements.  The  Fortran  unit  number  is  54.  The  records  are: 

1.  Record  1  contains  values  of: 


IEXPL,  I  DENS,  MM,  NA,  NB 
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Record  2  contains  values  of: 


VP1,  VS1,  R01 ,  VP 2,  VS2,  RQ2 


3.  Record  3  contains  all  values  of  array  VP32.  The  first 
dimension  of  tnis  array  is  the  value  of  MM.  The  second 
dimension  is  the  value  of 


N8NDY  =  NB  -  NA  +  2 
=  NB  -  NA  +  3 


when  ITRAN  =  1  or  ITRAN  =  2 
otherwise 


4.  Record  4  contains  all  values  of  array  VS32.  Dimensions  are 
the  same  as  VP32. 

5.  Record  5  is  used  only  when  IDENS  =  1.  It  contains  all  values 
of  array  R03.  The  dimensions  are  the  same  as  VP32. 

C.  Source  File 

This  file  is  created  by  FDSORS,  and  is  read  by  FINDIF.  It  is  a 
sequential,  unformatted  file,  created  with  binary  WRITE 
statements.  The  Fortran  unit  number  is  50.  The  records  are: 

1.  Record  1  contains  values  of: 

IEXPL,  I  FLAT,  IKELLY,  VP1,  KSTRT,  KK,  IRECT 

2.  Record  2  contains  values  of: 

ND,  MSW,  NSW,  DELT,  DELZ,  DELR,  NSORCE,  PLSWID,  TSWAVE 

3.  Record  3  contains  values  of  SORSA  for  the  first  time  step. 

The  dimension  of  SORSA  is: 

MM  if  IFLAT  =  1 

4*MSW  +  2*NSW  -  8  if  IFLAT  =  0 

4.  Record  4  contains  values  of  SORSB  for  the  first  time  step. 

The  dimension  is  the  same  as  for  SORSA. 

Records  3  and  4  are  repeated  for  each  time  step.  Total  number  of 
these  two  records  is  KK  -  KSTRT  +  1. 

D.  Snap-shot  Files 

Snap-shot  files  are  created  by  FINDIF.  They  are  sequential, 
unformatted  files,  created  with  binary  WRITE  statements.  The 
Fortran  unit  number  is  3.  The  first  snap-shot  file  is  created  at 
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time  step  KhAkK.  Thereafter,  a  snap-shot  file  is  created  every 
KMINC  timesteps,  until  the  final  timestep  (KK)  is  reached.  The 
names  of  these  files  are  created  by  the  program.  The  first  5 
characters  are  the  FILEID.  The  next  4  characters  are  the  numeric 
value  of  the  timestep.  The  file  extent  is  .SNS.  The  contents  of 
these  files  follow: 


1.  The  first  record  contains  values  of  all  integer  variables  in 
COMMON  block  MODPAR.  These  variables  are: 

LABEL(20) ,  IRECT,  IDENS,  I  FLAT,  IKELLY,  IEXPL,  IVERT,  ITRAN, 
ISORB,  NA,  NB,  ND,  NSW,  MSW,  MM,  NN,  KK,  KKSORS,  KOUTST, 
KOUTEN,  KINC,  NOUTST,  NOUTEN,  NINC,  MOUTST,  MOUTEN,  MINC, 
KSTRT,  KMARK,  KMINC,  NSORCE,  KCNT,  NSAkP 

2.  The  second  record  contains  values  of  all  real  variables  in 
COMMON  block  MODPAR.  These  variables  are: 


VP1,  VP2,  VPT,  VS1,  VS2 ,  VST, 
PLSWID,  TSWAVE,  TIME 


R01 ,  R02 ,  ROT, 


DELT,  DELR,  OELZ 


3.  The  third  and  last  record  contains  values  of  vertical 
displacement  over  the  entire  grid  (MM  *  NN  values). 


Time  Series  File 


This  file  is  created  by  Program  FINDIF.  It  is  a  sequential, 
unformatted  file,  created  with  binary  WRITE  statements.  The 
Fortran  unit  number  is  2.  The  name  of  the  file  is  the  same  as 
FILEID,  with  an  extent  of  .TST.  The  contents  of  the  time  series 
file  are  controlled  by  the  value  of  IVERT.  If  IVERT  =  0,  the  file 

will  contain  values  of  pressure;  if  IVERT  =  1,  it  will  contain 

values  of  vertical  displacement;  if  IVERT  =  2,  it  will  contain 

values  of  horizontal  displacement. 

The  number  of  individual  time  series  is  controlled  by  values  of 
MOUTST,  MOUTEN,  MINC,  and  NOUTST,  NOUTEN,  and  NINC.  One  series 
will  be  generated  for  every  MINC'th  range  grid  between  MOUTST  and 
MOUTEN,  and  every  NINC'th  depth  grid  between  NOUTST  and  NOUTEN. 
Each  time  series  will  include  every  KlNC'th  time  step  between 
KOUTST  and  KOUTEN. 

The  contents  of  the  records  in  the  file  are: 


1.  The  first  record  contains  values  of  all  integer  variables  in 
COMMON  block  MODPAR.  These  variables  are  listed  in  the 
snap-shot  file  description. 

2.  The  second  record  contains  values  of  all  real  variables  in 
COMMON  block  MODPAR,  listed  in  the  snap-shot  file  description 


2 
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3.  The  third  record  contains  the  following: 

LABEL (I) ,1=1,5  first  20  characters  of  label 

MLOC  range  grid  value  of  series 

NLOC  depth  grid  value  of  series 

KCNT  number  of  values  in  series 

NSAM3  number  of  samples  per  second  in  series 

RANGE  horizontal  distance 

TIME  time  of  first  value,  in  seconds 

ADW\X  maximum  absolute  value  in  series 

4.  This  record  contains  the  time  series  described  in  record  3. 
Record  types  3  and  4  are  repeated  for  as  many  time  series  as 
requested. 
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V.  The  subprograms 

A  number  of  subprograms  are  included.  These  are  listed,  with  brief 
descriptions,  in  Table  III.  It  will  be  noted  that  there  are  6  different 
TSTEP  routines.  The  one  used  depends  on  the  values  of  IRECT,  VS1,  and 
ITRAN,  according  to  the  table  below. 


IRECT 

VS1 

ITRAN 

TSTEP 

1 

*  0.0 

0 

1 

1 

t  0.0 

1 

1 

I 

*  0.0 

2 

1 

I 

-  0.0 

0 

1 

1 

-  0.0 

2 

3 

0 

*  0.0 

0 

2 

0 

t  0.0 

1 

2 

0 

*  0.0 

3 

5 

0 

A  0.0 

4 

6 

0 

a  0.0 

0 

2 

0 

•  -  0.0 

1 

4 

0 

=.  0.0 

2 

4 

0 

»  0.0 

3 

5 

0 

-  0.0 

4 

6 

Combinations  other  than  these  have  not  been  implemented. 
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Table  III 
Subprograms  used. 


NAME 

Used  by 

Variable 

COMMON 

needed 

Function 

LOGOUT 

all 

none 

Creates  first  part  of  Log  file 

RDMPAR 

all 

none 

Reads  model  parameter  file 

DATIM 

all 

none 

Access  current  date  and  time  from 
operating  system,  and  writes  them. 

OPNCOM 

FDPREP 

none 

Opens  files  for  COMMON  specifications 

BESJ 

FOSORS 

none 

Computes  Bessel  function  values  JO  or  J1 

BESY 

FDSORS 

none 

Computes  Bessel  function  values  YO  or  Y1 

COOL 

FDSORS 

none 

Does  FFT 

KELLY7 

FDSORS 

COMSOR 

Finds  displacement  due  to  Kelly  source 

SOURCE 

FDSORS 

none 

Finds  displacement  due  to  Alterman  source 

ABSORB 

FINDIF 

C0MFD8 

Does  absorbtion  calculation 

CONIT 

FIND1F 

C0MFD8 

Initializes  variables  and  arrays 

F  SET  UP 

FINDIF 

none 

Opens  files 

SCPOUT 

FINDIF 

none 

Outputs  snap-shot  file. 

TSOUT 

FINDIF 

none 

Outputs  time  series  file. 

TIMIT 

TIMSHO 

FINDIF 

none 

Initializes  system  timer 

Entry  to  TIMIT,  finds  elapsed  time. 

TSTEP1  ^ 

TSTEP2 

TSTEP3 

TSTEP4 

TSTEP5 

TSTEP6  J 

i  FINDIF 

C0ffD8 

These  routines  compute  the  displacements 
for  one  time-step. 
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VI.  Implementation  Notes 

The  programs  have  been  written  in  Fortran  77  for  a  VAX  11/780  under 
VMS.  An  effort  has  been  made  to  avoid  non-standard  features.  Specifically, 
all  identifiers  are  6  characters  or  less,  and  the  DO  WHILE  feature  is  not 
used.  The  programs  all  use  the  BLOCK  IF  capability,  and  do  use  character 
variables.  If  the  INCLUDE  statement  is  not  available,  programs  and  some 
subroutines  will  have  to  be  edited  to  change  some  of  the  array  dimensions. 

The  various  TSTEP  routines  are  rather  long,  so  we  have  created  a  subroutine 
library  which  includes  dummy  versions  of  them.  When  we  load  FINDIF,  we 
specifically  include  the  TSTEP  needed,  and  satisfy  the  unused  references  to 
the  other  TSTEP  routines  from  the  library.  It  may  be  necessary  at  some 
installations  to  use  some  other  technique. 

An  effort  has  been  made  to  isolate  all  input  and  output,  and  all 
non-standard  features,  in  subroutines.  The  following  routines  might  need 
changes: 

TIMIT  to  initialize  the  system  timer 

TIMSHO,  an  ENTRY  to  TIMIT,  finds  elapsed  CPU  time 

FSETUP  opens  output  Log  file,  source  file,  and  boundary  velocity  file 

SNPOUT  creates  name  of  snap-shot  file,  opens  file,  writes  data  in  file,  and 
closes  file 

TSOUT  creates  name  of  time  series  file,  opens  file,  writes  data  in  file,  and 
closes  file 

• 

RDMPAR  used  by  all  programs,  reads  the  model  parameter  file.  It  uses 
list-directed  READ  statements. 

DATIM  used  by  all  programs,  gets  current  date  and  time  from  operating  system. 
OPNCOM  used  by  FDPREP,  opens  files  for  new  COMMON  specifications. 
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VII.  Usage  on  the  VAX 

A.  General 

This  section  describes  how  the  system  is  set  up  to  run  on  the  VAX 
here  at  WHOI.  Changes  might  be  needed  for  other  computers.  In 
general,  ASSIGN  commands  are  required  for  input  files,  and  output 
file  names  are  created  by  the  programs  using  the  FILE  I D  from  the 
model  parameter  file  with  a  different  file  extent  (or  file  type) 
for  each  output  file. 

B.  Model  Parameter  File 

The  first  step  is  to  create  a  model  parameter  file.  The  contents 
of  the  records  are  listed  in  Section  IV  A,  and  descriptions  of  the 
variables  are  in  Table  I.  All  records  except  the  first  two  are  in 
free-field  format,  so  spacing  is  not  important.  Let  us,  in  this 
example,  call  this  file  MODEL. PAR,  and  let  us  use  M0DL1  as  the 
file  identifier. 

C.  Run  F DP REP 

Next,  Program  FDPREP  should  be  run.  This  program  does  the 
following: 

1.  Checks  model  parameters  for  consistency 

2.  Creates  two  COMMON  block  files: 

COM7 D8. FOR  used  by  FDBNY  and  FINDIF 
COMSOR .FOR  used  by  FDSORS. 

3.  Creates  a  Log  file  summarizing  its  activities.  In  this  case, 
the  file  will  be  M0DL1.LG1 

Before  running  FDPREP,  you  must  use  an  ASSIGN  command  to  connect 
the  model  parameter  file  to  unit  55.  Commands  to  run  FDPREP  are 

SASSIGN  MODEL. PAR  F0R055 

$RUN  FDPREP 

D.  Run  FDBNY 

Now  Program  FDBNY  can  be  compiled  and  run.  Here  again,  the  only 
ASSIGN  command  required  is  for  the  model  parameter  file.  The 
program  reads  this  file  and  creates  the  following  two  files: 


I 
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1.  The  file  defining  the  velocity  structure  in  the  transition 
zone.  The  name  of  this  file  will  be  MOOLI. BNY 

2.  A  log  file,  named  M0DL1.LG2 

The  commands  to  compile  and  run  FDBNY  are: 

SFORTRAN  FDBNY 
SLINK  FDBNY .FOLIB/LIB 
SASSIGN  MODEL. PAR  F0R055 
SRUN  FDBNY 

Run  FDSORS 

He  can  now  compile  and  run  Program  FDSORS.  This  Is  similar  to  the 
procedure  for  FDBNY.  The  files  created  are: 

1.  A  file  defining  the  energy  source  for  each  time  step. 

The  name  of  this  file  will  be  MODL1.SOR 

2.  A  log  file,  named  M0DL1.LG3 

The  commands  to  compile  and  run  FDSORS  are: 

SFORTRAN  FDSORS 
SLINK  FDSORS, FDLIB/LIB 
SASSIGN  MODEL. PAR  F0R055 
SRUN  FDSORS 


Subroutine  KELLY7  is  in  the  same  file  as  FDSORS. 

Run  FINOIF 

Finally,  we  can  compile  and  run  FINDIF.  He  must  also  compile  the 
version  of  TSTEP  which  will  be  needed.  This  information  Is 
Included  in  the  log  file  created  by  FDPREP.  The  commands  to  run 
FINDIF  are: 

SFORTRAN  FINDIF 

SFORTRAN  TSTEPn  use  the  selected  version 

SFORTRAN  FINSUB  this  file  Includes  ABSORB  and  CONIT 

SLINK  FINDIF, TSTEPn .FINSUB,  FDLIB/LIB 

SASSIGN  MODEL. PAR  F0R055 

SASSIGN  M0DL1.BNY  F0R054 

SASSIGN  M0DL1.S0R  F0R050 

SRUN  FINDIF 
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The  files  created  by  FINDIF  are: 

1.  Log  file,  named  M0DL1.LG4 

2.  A  time  series  file,  named  MODLl.TST 

3.  One  or  more  snap-shot  files.  The  names  includes  the  time 
step.  For  example,  a  snap-shot  file  of  time  step  1000  will 
be  MOOLI 1 000. SNS 

G.  Boundary  and  Source  Files 

It  is  not  necessary  to  create  new  boundary  and  source  files  for 
each  run  of  FINDIF.  However,  there  are  some  parameters  which  must 
be  the  same  as  those  used  to  create  the  files.  For  FDBNY,  these 
parameters  are  MM,  NA,  NB,  VP1,  VS1,  R01,  VP2,  VS2,  and  R02.  For 
FOSORS,  they  are  IFLAT,  IKELLY,  VP1,  KSTRT,  KK,  ND,  MSW,  IRECT, 
DELT,  OELZ,  DELR,  NSORCE,  PLSWID,  and  TSWAVE.  If  any  of  these 
parameters  in  FINDIF  is  not  the  same  as  that  used  to  create  the 
boundary  and/or  source  files,  an  error  message  is  output,  and 
FINDIF  terminates. 
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VIII.  Error  and  Diagnostic  Messages 

The  following  error  messages  are  included  so  the  conditions  causing  the 
errors  can  be  avoided. 

A.  FOPREP 

1.  The  first  two  messages  are  the  result  of  inconsistencies  in 
the  time  step  parameters.  KSTRT  and  KK  are  the  first  and  last 
time  steps.  KOUTST  and  KM4RK  are  the  first  time  step  for  time 
series  and  snap-shot  files  respectively.  KOUTEN  is  the  time 
step  for  the  end  of  the  time  series.  The  messages  are: 

KOUTST  = _ or  KMARK  = _  IS  LESS  THAN  KSTRT  = _ _ 

KOUTEN  = _  IS  GREATER  THAN  KK  = _ 

2.  When  IFLAT  =  0,  NSW  must  be  an  odd  number.  If  not,  the 
following  message  will  be  output 

NSW  IS  NOT  AN  000  NUMBER  -  JOB  STOPPED 

3.  Invalid  combinations  of  parameters  can  produce  the  following: 

INVALID  COMBINATION:  ISORB  » _  ITRAN  = _ 

INVALID  COMBINATION:  IRECT  = _ VS1  = _ ITRAN  = _ 

4.  The  last  is  a  warning.  The  program  will  continue. 

WARNING  -  ABSORBING  BOUNDARIES  NOT  CODED  FOR  THIS  COMBINATION. 

B.  FINDIF 

1.  The  program  will  not  run  if  the  parameters  used  to  create  the 
boundary  and/or  source  files  do  not  match  the  current 
parameters.  The  messages  are: 

BOUNDARY  PARAMETERS  DO  NOT  MATCH 

SOURCE  PARAMETERS  DO  NOT  MATCH 

In  each  case,  the  offending  parameter  values  are  output. 

2.  The  next  message  can  appear  only  when  the  model  parameter  file 
is  changed  after  running  FDPREP,  or  the  user  does  something 
stupid. 


IWALID  COMBINATION  OF  PARAMETERS:  IRECT 


ITRAN 


VS1 


File  errors: 


ERROR _  OPENING  FILE _ 

SNAPSHOT  FILE  ERROR.  CODE  = _  KMARK  = _ 

ERROR  ON  TIME  SERIES  FILE.  CODE  = _ 

As  can  be  seen,  in  each  case  the  error  code  is  included  in  the 
message. 
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IX.  Sample  Runs 

A  group  of  samolo  runs  was  made,  using  the  same  basic  parameters;  the 
only  changes  in  the  parameter  file  were  those  required  by  the  use  of  different 
options.  The  contents  of  the  first  model  parameter  file  are  listed  with  the 
output  of  the  first  run.  The  changes  to  create  the  other  parameter  files  are 
given  in  Table  IV.  The  sample  run  of  the  first  model  follows,  step  by  step, 
the  instructions  in  Section  VII.  Simple  contour  plots  of  the  snap-shot  files 
and  seme  time  series  plots  are  included. 

It  will  be  noticed  that  the  snap-shots  for  STEV8  and  STEVD  {cylindrical 
coordinates  and  a  square  source)  include  some  anomalies  in  the  source  area. 
This  noise  is  a  result  of  the  small  size  of  the  model  and  is  not  present  in 
larger  models. 

A.  Model  Parameter  File 

Below  is  a  listing  of  the  model  parameter  file  used  in  the  first 
run.  As  can  be  seen,  the  name  of  the  file  is  STEV1.PAR,  and  the 
FILEID  is  also  STEV1.  These  two  names  do  not  have  to  be  the  same, 
although  it  saves  some  confusion  if  they  are. 


STEV1  FINAL  TEST  FOR  RELEASE 
STEOl 

It  It  It  1?  It  It  Ot  1 

72 1  72 1  300/  1/  0,003/  0,04/  0.04 
It  300/  It  72/  37»  37 
2 1  5 1  1/  100»  100 

1.5/  0.0/  1.0/  2.5/  1.4/  1.5/  17/  47/  0,0/  0.0/  0.0 
-11 /  0/  0 
i  n  -  o  .  ft  - 

A  V  /  W  V  ♦  M  /  W  «  W  t»  /  t. 
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Table  IV 

Parameter  Values  for  Test  Runs 


Parameter 


f  i  le 

TSTEP 

STEV1 

STEV2 

STEV3 

STEV4 

STEV5 

STEVA 

STEVB 

STEVC 

STEVD 

routine 

1 

1 

2 

3 

4 

1 

2 

3 

4 

IFLAT 

1 

1 

1 

1 

1 

0 

0 

0 

0 

NO 

-11 

-11 

-11 

-11 

-11 

10 

10 

10 

10 

MSW,NSW 

0,0 

0,0 

0,0 

0,0 

0,0 

5,11 

5,11 

5,11 

5,11 

ITRAN 

0 

0 

0 

2 

2 

0 

0 

2 

2 

NA,NB 

17,47 

17,47 

17,47 

30,30 

30,30 

17,47 

17,47 

30,30 

30,30 

IRECT 

1 

1 

0 

1 

0 

1 

0 

1 

0 

ISORB 

1 

0 

0 

1 

0 

1 

0 

0 

0 

Boundary 

file 

Source 

STEV1 

STEV1 

STEV1 

STEV4 

STEV4 

STEV1 

STEV1 

STEV4 

STEV4 

file 

STEV1 

STEV1 

STEV3 

STEV1 

STEV3 

STEVA 

STEVB 

STEVA 

STEVB 

FINDIF 

CPU  time 
(Minutes) 

5:50 

5:41 

6:38 

2:52 

3:43 

6:15 

6:35 

2:46 

3:42 
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S.  Run  Program  FDPREP 

A  listing  of  the  log  file  is  included 


i  ASS  I  Gil  STEVl.FAR  F0R055 
»  HUN  f  DTRCP 

FDPREP  COMPLETED  — NO  ERRORS 
i  IYFE  3  T  E  V 1  .LGi 


PROGRAM  FDPREP  VERSIOH  OF  20 -SEP-83 

DATE  AMD  TIME  OF  RUM  23-SEP-S3  13150 

FILE  ID  IS  5 TEL1 1 
VALUES  OF  INPUT  PARAMETERS: 


IRcCT 

I  TRAN 

=  1 

=  0 

I  DENS  =  1 

ISORB  =  1 

I  FLAT 

=  1  IKELLY  = 

1  IEXPL  = 

1  I  VERT 

MM  = 

92 

NN  =  72 

KK  ■= 

300  K  S  T  F;  T  - 

. 

Hh  = 

17 

MB  =  4  7 

MD  = 

-11  MSN  = 

0  NSW  = 

0 

DELT  - 

0. 

008 

DELR  ~ 

0.0  40 

DELZ  = 

0.0  40 

OUTPUT  SPAMS 

TIMES 

AND  I 

1 

ilCREHENTS: 

300 

2 

RANGES 

1 

72 

o 

DEPTHS 

37 

37 

1 

SNAPS 

100 

100 

VELOCITIES! 

VP1  =  1.500 

VS1 

0.000 

E01  = 

1.000 

VP  2  =  2.500 

V32 

=  1.400 

U02  = 

l  .500 

V  P  T  -  0.000 

VST 

-  0 . 000 

ROT  - 

0.000 

FELLY  SOURCE  PARAMETERS  t 

i!  GORGE  -  10  Pl.SUID  =  50.0000  TSUAVE  =  0.3272 


FIMDIF  WILL  USE  SUBROUTINE  TSTEPl 

N  A.  INTERCEPTED  A  T  TIME  STEP  93 

Mil  INTERCEPTED  NO  EARLIER  THAN  TIME  STEF",  1S1 

HM  INTERCEPTED  HO  EARLIER  1  HAM  TIME  STEP:  101 


DIMENSION  FOR  ARRAYS 


r  u 

A  AND 

EC  7  2 

72  2 

r  OP 

AR  AND 

br  :  i 

l  2 

rojv 

f'OF: 

C»  D1 
an:  < 

(  72*2) 

72,2;  ru.C 

<  7  7 ) 

rc.F, 

T.  I  Gil  i  n 

*  r  I  til .  1  i  -  ,  i‘i 

A  M  1  ,  R  1 5  '<  ii  i : 

l  l .  i  \' 

r.o  r .  " 

nT?..r:i)T.,:;i , 

r.orz.u: 

i  UK 

aila  i  ii : 

1  5 

:  iso 

r  Hu 

C  U  SA 

AMU  CURSE'. 

/ 

r  *.  . 

M  |  •  ?  “j 

MD  VS 32: 

)  2  ^  Z 

1  43 
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C.  Compile,  link  and  run  FDBNY 

Since  this  was  done  immediately  after  runnino  FDPREP,  we  did  not 
need  to  repeat  the  ASSIGN  command  for  STEV1.PAR.  A  partial  listing 
of  the  log  file  is  included. 


i  T GRTRAN  FDBNY 
t  LINK  EBONY  fFDL IB/LIB 
l  RUN  FDBNY 
FORTRAN  STOP 
1  TYPE  STFV 1  . LG2 


PROGRAM  FDBNY  VERSION  OF  9 -JUL-33 

DATE  AMD  TIME  OF  RUN  23-SEP-33  13151:51 

FILE  ID  IS  STEV1 
VALUES  OF  INPUT  PARAMETERS*. 


IRSCT 

I  TRAN 

-  1 

-  0 

IDENS 

ISORB 

-  1 

=  1 

IFLAT 

=  1  IKELLY  = 

1 

IEXF'L  =  1  IVERT 

MM  = 

72 

NN  = 

72 

KK  = 

300  KSTRT  = 

1 

NA  = 
BELT  = 

17  NB  = 

0.008 

47 

ND  = 
DELR  = 

-11  MSU  = 

0.040 

0 

MSU  =  0 

DELZ  =  0.040 

VELOCITIES*. 

VP1  =  1.500  VS1  = 

VP 2  •••  2.500  VS2  = 

VFT  =  0.000  VST  = 

VP * *2  IN  TRANSITION  Zi 

0.000 
1.400 
0.000 
ONE »  (It 

R01  = 
R02  = 
ROT  = 
J)  ,  I  = 

1.000 

1 . 500 
0.000 
If  72 

r  3  J 

*  1. 

33  • 

2.230 

2.250 

2.250 

2.250 

2,250 

2.250 

2.250 

2 . 250 

2.250 

2.441 

2.441 

2.441 

2.441 

2.441 

2.441 

2,441 

2 . 4  A  1 

2.441 

2.510 

2.510 

2.540 

2.340 

2.540 

2.540 

2.5  40 

2.340 

2.540 

2.641 

2.641 

2  .  -i  1 1 

2.441 

2.641 

2.641 

2.641 

2.641 

2.641 

2.7  43 

2.743 

2.743 

2.743 

2.743 

2,743 

2.743 

2.743 

2.743 

2.343 

2.343 

2.343 

2.313 

2.318 

2.843 

2.343 

2.843 

2.843 

2.951 

2.754 

2.954 

2.954 

2.954 

2.954 

2.954 

2.954 

2.954 

3.063 

o ,  0  6  o 

3.063 

3  ♦  0  o  o 

3  •  063 

^  O 

3 . 063 

3.063 

3.063 

3.173 

3.173 

3.173 

3.173 

3,173 

3.173 

3.173 

3.173 

3.173 

3.235 

3.285 

3.235 

3.285 

3.233 

o  .  205 

3.235 

3.235 

3.285 

3 . 399 

o  .  39  9 

3.377 

3.399 

3.399 

3.399 

3.399 

3.399 

3.399 

3.516 

3.516 

3.516 

3  .  o  1  o 

3.316 

3,516 

o  .  -j  1 6 

3.516 

3.514 

3.634 

3.634 

3.03  1 

8  •  634 

3.634 

3 .634 

3.634 

3.634 

3.634 

3.751 

3.754 

3.754 

3  .  >’54 

3.754 

3.754 

3.754 

3.754 

3.754 

3.376 

3.376 

3.376 

3.376 

3.376 

3.876 

3. 376 

3.376 

3.376 

4.000 

■1.000 

4.000 

•1.000 

4 .000 

■1.000- 

1.000 

4.000 

4.000 

4.126 

4.126 

4.126 

1.126 

4.126 

1  .  126 

4.126 

1.126 

4 .126 

1.254 

1.251 

4.231 

4.234 

1.251 

1.254 

•1  .  254 

4.25  4 

'4.254 

4.384 

4.304 

4.331 

4 . 304 

1.334 

4.334 

'..33 1 

•1.504 

A  .  33  l 

4.516 

4.516 

4.516 

4.516 

4.516 

4.516 

4.516 

4.516 

1.516 

4 . 417 

•1.447 

4 . 44? 

4.649 

4.649 

4.649 

4.649 

1.6  17 

•1.649 

4.735 

4  .  785 

4 . 735 

4.735 

4.735 

4.735 

4.735 

4.785 

4 . 785 

4.927 

.»  O  IT 

I  *  •  s  uJ 

1 .723 

4.723 

4.723 

1 . 723 

1.727 

4.723 

4.923 

3.063 

5 .063 

5.063 

5  .  0o3 

5.063 

5  ■  0  0  3 

3.06  3 

5.063 

3.0*3 

5 . 20  4 

5.201 

5.201 

5 .20  4 

5 . 2  4 

*2.204 

5 . 2  0  4 

5.204 

5 .  204 

5  .  o  4  •  l 

5.348 

5.318 

5.318 

5.343 

D.  Compile,  link,  and  run  FDSORS 


Again,  we  have  included  a  listing  of  the  log  file. 


$  FORTRAN  FDSORS 
f  LINK  FDSORS fFDLIB/L IB 

*  PUN  FDSORS 
FORTRAN  STOP 

*  TYPE  STEV1.LG3 


PROGRAM  FDSORS  VERSION  OF  12-SEP-33 

DATE  AND  TIME  OF  RUN  23-3EP-S3  13:5 

FILE  ID  IS  3TEV1 
VALUES  OF  INPUT  PARAMETERS} 


IRECT 

ITRAN 

•-  1 

=  0 

I  DENS  =  1 

ISORB  *  1 

IFLAT 

=  1  IK'ELLY  = 

1 

IEXPL  =  1  IVERT 

MM  = 

72 

NN  =  72 

KK  = 

300  KSTRT  = 

1 

NA  ~ 

17 

NB  =  47 

ND  = 

-11  MSW  = 

0 

NSW  =  •  0 

BELT  = 

0. 

003 

DELR  = 

0.040 

DELZ  =  0.040 

OUTPUT  SPANS  AND  INCREMENTS} 


TIMES 

1 

300 

2 

RANGES 

1 

72 

5 

DEPTHS 

37 

37 

1 

SNAPS 

100 

100 

velocities: 

ypi  = 

1.500 

VS1  = 

0.000 

R01  = 

1.000 

VP2  = 

2.500 

VS2  = 

1  .  400 

RQ2  - 

1.500 

VPT  -= 

0.000 

VST  = 

0.000 

ROT  = 

0.000 

KELLY  SOURCE  PARAMETERS! 


NSDRCE  =  10  PLSUID  =  50.0000  T3UAVE  =  0.5292 
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E.  Run  FINDIF 

Finally,  we  compile,  link,  and  run  FINDIF.  We  use  TSTEP1  because 
Program  FDPREP  informed  us  which  TSTEP  routine  would  be  needed.  As 
can  be  seen.  Program  FINDIF  displays  the  elapsed  time  and  CPU  time 
every  50  time  steps.  The  log  file  is  included  on  the  next  page. 


$  FORTRAN  FINDIF 
t  FORTRAN  TSTEP1 
i  FORTRAN  FINSUB 
»  ASSIGN  STEV1.SOR  F0R050 

*  ASSIGN  STEV1.BNY  F0R054 

*  LINK  FINDIF  »TSTEP1 »FINSUB»FDLIB/LIB 

*  RUN  FINDIF 


PERFORMANCE  STATISTICS 

AT  TIME  STEP 

50 

elapsed:  00:01:03.53 

cpu:  0:01:01.03 

bufio: 

n 

dirio:  3 

faults:  59 

PERFORMANCE  STATISTICS 

AT  TIME  STEP 

100 

elapsed:  00:02:13.00 

CPU!  0:02:03.30 

bufio: 

3 

dirio:  10 

faults: 

10? 

PERFORMANCE  STATISTICS 

AT  TIME  STEP 

150 

elapsed:  00:03:27.43 

cpu:  0:03:05.32 

bufio: 

11 

dirio:  13 

faults: 

261 

PERFORMANCE  STATISTICS 

AT  TIME  STEP 

200 

ELAPSED:  00: 04:40.84 

cpu:  0:01:07.5? 

bufio: 

17 

dirio:  20 

faults: 

230 

PERFORMANCE  STATISTICS 

AT  TIME  STEP 

250 

elapsed:  00: 06: 20. 06 

cpu:  0:05:10.19 

bufio: 

20 

dirio:  23 

faults: 

283 

PERFORMANCE  STATISTICS  AT  TIME  STEP  300 

elapsed:  00:07:52.30  cpu:  o:os:  12.73  dufio:  26  dirio:  29  faults:  285 

FORTRAN  STOP 
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t  TYPE  S  TCV  l  .1  G4 


PROGRAM  FIMDih  vtfxalun  Or  22  ii w G  u3 

CATC  AMD  TIME  OF  RUN  23-SEP-83  13 159  *.56 

FILE  ID  IS  STEV1 
VALUES  OF  INPUT  FAKAMEILKSi 


IRECT  a  1 

IDENS  =  1  IFLAT  =  1 

I KELLY 

1 

IEXFL  r  i  I VERT 

I  TRAN  =  0 

ISORD  -  1 

MM  =  72 

NN  =  72  KK 

a  300 

KSTRT  a 

1 

NA  =  17 

NB  =  47  ND 

=-  -11 

HSU  = 

0 

NSU  =  0 

CELT  =  0.008 

DELR  =  0. 

040 

IiELZ  =  0.040 

OUTPUT  SPANS 

AND  increments: 

TIMES 

1  300 

•i 

RANGES 

1  72 

s 

DEPTHS 

37  37 

1 

SNAPS 

100 

100 

velocities: 

VP1  =  1.500 

VS1  =  0.000 

R01  - 

1.000 

VP2  =  2.500 

VS2  =  1.400 

R02  = 

1.500 

VFT  =  0.000 

VST  =  0.000 

ROT  = 

0.000 

KELLY  SOURCE 

parameters: 

NSORCE  =  10 

PLSWID  = 

50.0000  TSUAVE  = 

0.5292 

2.250 

2.441 

2.540 

2.641 

2,743 

2.848 

• 

2.954 

3.043 

3.173 

3.285 

3.399 

3.516 

a  * 

3 . 754 

♦  S  /  4 

4.000 

4 . 126 

4.254 

4.334 

4,514 

4.64? 

4.785 

4.923 

5.063 

5.20  4 

5.348 

5 . 493 

5.641 

5.790 

5.941 

<4.095 

4.250 

6.250 

0.000  • 

0.008 

0.017 

0.031 

0.048 

0.069 

0.094 

0.122 

0.155 

0.191 

0.232 

0.276 

0.323 

0.375 

0.431 

0.490 

0.553 

0.620 

0.691 

0.766 

0.344 

0.926 

1.013 

1 . 102 

1.194 

1.294 

1.395 

1.501 

1.610 

1.723 

1 .  .33? 

1.940 

1.940 

1.000 

1.031 

1.047 

1.063 

1.078 

1 .094 

1  .  10? 

1.125 

1.141 

1.156 

1.172 

1.188 

1.203 

1.219 

1.234 

1 .250 

1.266 

1.231 

1.297 

1.313 

1.328 

1.344 

1.357 

1.375 

1.391 

1 . 406 

1 .  422 

1.433 

1.453 

1.469 

1.434 

1.500 

1.500 

VP  3c  *2  til 

TRANS  I 

7 ION  ZONE/  <Ir 

J)  r  I-=l  , 

72 

8 

J^l 

i 

33 

2.250 

2 . 250 

2.250 

2.250 

2.250 

2.250 

2 

,250 

fc  » 

250 

<.  «  b  jO 

2.441 

2 . 4  4  l 

2.441 

2.441 

2.441 

2.441 

2.441 

441 

1  51 

2  .  2  ■  '*  c* 

2.540 

2.540 

,  j 40 

2.340 

2.540 

2.54  0 

2.540 

540 

2  . 

oil 

2  *  i  U 

1.641 

?.4U 

2.611 

2.441 

2.641 

2 . 4  4  i 

2.611 

n 

7  13 

Z  • 

2  15 

%  *7  1  *T 

2.745 

2 ,743 

2.743 

2.743 

2.713 

2 . 7  'i  3 

.  n  i » 

7  15 

6  1 5 

"4  *  •  "' 

ft-  »  \J  i  v. 

2 .542 

’2 . 345 

2.340 

2.340 

2.340 

2.954 

2.954 

■> 

954 

n 

/  5  4 

?.  ,-za 

2.754 
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The  following  pages  contain  contour  plots  of  the  snap-shot  files.  A 
few  words  of  explanation  about  the  format  of  these  plots  will  help  in 
understanding  them.  The  FI  LEI D  and  time  step  number  are  in  the  lower  right 
corner  of  each  plot.  The  number  in  the  lower  left  corner  is  the  scaling 
factor  used  for  the  plot,  which  is  selected  by  the  plot  program.  In  some 
cases,  we  have  used  the  same  scaling  for  all  timesteps,  and  in  others,  we  have 
used  a  different  scale  for  each  timestep.  The  dashed  lines  indicate  negative 
contours. 
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X.  Notes  on  Parameter  Specifications 

The  values  used  in  the  parameter  file  cannot  be  chosen  arbitrarily. 
For  stability  of  the  code 


at  < 


minUr.AZ 


.  jnn 

p 


where  a  and  s  are  congressional  and  shear  wave  speeds.  For  acceptable  grid 
dispersion  one  should  choose  at  least  ten  grid  points  per  wavelength  at  the 
upper  half  power  frequency.  That  is 


AX 


,  a1/2 
^  ~nr 


For  the  Kelly  source  (the  Alterman  source  is  not  recommended)  the  center 
frequency  for  displacement  is  f0  =  0.318^£  and  the  upper  half  power 
frequency  is  f+1/2  =  0.458\/f  . 

Codes  without  density  variations  have  not  yet  been  written  so  IDENS  is 
always  one.  Also  the  absorbing  boundaries  only  work  for  rectangular 
co-ordinates.  This  program  contains  only  explicit  codes  so  IEXPl  is  always 
one.  When  I FLAT  equals  one  the  source  is  introduced  along  the  top  of  the  grid 
as  if  the  source  were  above  the  box  (ND  is  negative)  and  NSW  and  MSW  are  not 
used. 


When  IFLAT  equals  zero  the  source  is  introduced  on  the  edges  of  an  MSW 
X  NSW  box  at  the  axis  of  symmetry  centered  at  depth  ND  (positive).  TSWAVE 
should  be  chosen  so  that  the  Initial  disturbance  at  the  grid  is  small  (~  10"°). 
Examples  of  a  number  of  different  options  are  included  in  the  test  runs. 


i 
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APPENDIX 

Performance  Evaluation 


The  calculation  of  synthetic  seismograms  by  finite  difference  methods 
requires  significant  amounts  of  CPU  time  and  computer  memory.  Throughout  the 
program  development  we  have  tried  to  minimize  the  resources  required,  to 
enable  us  to  run  larger  and  more  complex  models  at  reasonable  cost.  During 
March,  1981  a  study  was  conducted  to  determine  how  the  relation  between 
working  set  size  and  program  size  affects  CPU  time  and  page  faults.  Three 
different  sized  models  (1.2  Mbyte,  1.8  Mbyte,  and  2.7  Mbyte)  were  run  with 
working  set  sizes  varying  from  400  to  4000  pages.  Table  V  summarizes  the 
resul ts. 

It  will  be  noted  that,  for  each  size  model,  there  is  a  critical 
working-set  size,  where  the  number  of  page  faults  is  drastically  reduced. 
Nothinq  is  gained  by  increasing  the  working-set  size  above  this  critical 
size.  It  will  also  be  noted  that  the  difference  between  the  program  size  and 
the  critical  working-set  size  is  the  same  for  each  model.  This  constant,  we 
can  assume,  is  the  size  of  the  code  which  is  only  executed  once,  while  the 
critical  working-set  size  includes  the  arrays  used  in  the  computation  and  the 
repeated  code.  Since  the  largest  model  in  these  tests  is  still  a  small  model, 
and  the  4000  page  working-set  required  was  pushing  the  limit  of  the  VAX,  it  is 
obvious  that  some  drastic  steps  would  be  required  to  run  even  moderately  large 
models  economically  on  the  VAX. 

A  second  group  of  tests  was  made  to  determine  the  effect  of  using 
larger  arrays  than  needed  for  the  model.  A  model  which  would  have  required 
1.5  Mbytes  was  run  in  a  program  of  1.8  Mbytes  and  one  of  2.7  Mbytes.  The 
results  are  given  in  Table  VI.  Although  the  working  sets  used  were  below  the 
critical  size  in  each  case,  it  can  be  seen  that  CPU  time  was  reduced 
substantially  by  using  the  smaller  arrays.  This  is  the  reason  for  recompiling 
and  relinking  the  programs  for  each  model. 


Relation  between  program  size. 

working- set  size. 

and  page  faults. 

PEAK 

DEFAULT 

PAGE 

WORKING 

WORKING 

FAULTS 

CPU  TIME 

SET  SIZE 

SET  SIZE 

2997 

6.92 

400 

400 

2258 

6.54 

600 

600 

2229 

6.55 

800 

800 

1271 

6.08 

1000 

1000 

1289 

6.03 

1129 

1200 

1240 

6.07 

1130 

1400 

1238 

6.17 

1131 

1600 

1227 

7.13 

1141 

1800 

1353 

6.99 

1141 

2000 

1.2  Mbyte  model 

(2392  pages) 

PEAK 

DEFAULT 

PAGE 

WORKING 

WORKING 

FAULTS 

CPU  TIME 

SET  SIZE 

SET  SIZE 

33705 

75.81 

500 

500 

28458 

74.53 

1000 

1000 

25887 

72.97 

1500 

1500 

25790 

73.06 

2000 

2000 

2443 

63.24 

2318 

2500 

2442 

63.63 

2322 

3000 

2442 

63.63 

2326 

3500 

2442 

64.71 

2331 

4000 

1.8  Mbyte  model 

(3595  pages) 

PEAK 

DEFAULT 

PAGE 

WORKING 

WORKING 

FAULTS 

CPU  TIME 

SET  SIZE 

SET  SIZE 

41637 

134.81 

500 

500 

41121 

135.45 

1000 

1000 

36857 

134.73 

1500 

1500 

30381 

135.60 

2000 

2000 

30695 

133.77 

2500 

2500 

34092 

130.09 

3000 

3000 

41383 

124.15 

3500 

3500 

4189 

102.17 

3998 

4000 

2.7  Mbyte  model  (5267  pages) 


Table  VI 


PEAK 

DEFAULT 

PAGE 

WORKING 

WORKING 

FAULTS 

CPU  TIME 

SET  SIZE 

SET  SIZE 

35537 

83.91 

200 

200 

30114 

84.36 

400 

400 

31200 

81.57 

600 

600 

29764 

82.33 

800 

800 

28103 

84.62 

1000 

1000 

25655 

85.16 

1200 

1200 

23946 

84.76 

1400 

1400 

22867 

83.98 

1600 

1600 

21394 

82.61 

1800 

1800 

23176 

83.68 

2000 

2000 

1.5  Mbyte  model  in 

2.7  Mbyte  program  (5267  pages) 

PEAK 

DEFAULT 

PAGE 

WORKING 

WORKING 

FAULTS 

CPU  TIME 

SET  SIZE 

SET  SIZE 

35660 

78.57 

200 

200 

34002 

78.41 

400 

400 

33617 

77.99 

600 

600 

31749 

76.35 

800 

800 

26455 

74.89 

1000 

1000 

25930 

75.20 

1200 

1200 

25590 

73.97 

1400 

1400 

25913 

74.66 

1600 

1600 

25666 

74.19 

1800 

1800 

25767 

73.84 

2000 

2000 

1.5  Mbyte  model  In  1.8  Mbyte  program  (3595  pages) 
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As  a  result  of  these  studies  It  became  clear  that  more  computing  power 
would  be  necessary  to  run  "real -world"  models.  The  first  step  was  to  increase 
the  physical  memory  of  the  VAX  from  2.5  to  4.0  megabytes.  The  performance  when 
running  small  to  moderate  models  was  improved,  but  it  was  still  not  reasonable 
to  run  large  models. 

Array  Processors  were  then  explored  as  a  means  of  improving  performance. 
The  Floating  Point  Systems  family  of  array  processors  was  chosen  as  a  possibil¬ 
ity  and  a  series  of  bench-marks  was  run.  Table  VII  summarizes  the  results.  It 
is  obvious  that  the  VAX/AP  performance  is  degraded  by  the  excessive  amount  of 
I/O  between  the  VAX  and  the  A P.  In  fact,  even  moderate-sized  models  are  I/O 
bound.  Until  the  I/O  bandwidth  is  increased,  or  sufficient  AP  memory  is  avail¬ 
able,  the  array  processor  is  not  the  solution. 

The  only  remaining  option  seemed  to  be  the  use  of  a  super  computer.  In 
December  1982,  a  contract  was  signed  with  Control  Data  Corporation  to  estab¬ 
lish  a  CYBERNET  SERVICES  link  to  a  CYBER-205.  The  initial  code  conversion  was 
performed  by  Stu  Gould  of  C.D.C.;  he  also  provided  vectorized  code  for  two 
main  loops  in  the  TSTEP  routine.  Work  began  in  February  1983  in  Woods  Hole 
via  a  300  BAUD  dial-up  phone  line.  After  operational  details  had  been  ironed 
out,  successful  runs  of  both  the  scalar  and  vector  codes  were  made.  A  com¬ 
parison  of  CYBER-205  and  VAX  timing  and  cost  is  shown  in  Table  VIII. 

It  must  be  realized  that  while  these  figures  provide  a  "ball-park"  idea 
of  efficiency,  they  are  not  exactly  comparable.  VAX  charges  reflect  only  CPU 
cost  while  CYBER  charges  are  total  costs.  This  error  is  about  10  .  The  vec¬ 
tor  code  becomes  more  efficient  as  the  model  size  increases,  leading  us  to 
expect  more  favorable  costs  for  a  larger  model.  Finally,  the  CYBER-205  runs 
were  done  during  prime  shift,  which  adversely  affects  wall -clock  time,  while 
the  VAX  run  was  non-prime  shift  on  a  dedicated  machine,  optimizing  wall-clock 
time. 


Since  that  time,  work  has  been  progressing  toward  implementing  a  later 
version  of  the  programs  in  the  CYBER-205.  Preliminary  tests  Indicate  that  the 
unvectorized  code  runs  about  60  times  faster  on  the  CYBER  for  large  models. 

We  hope  for  even  further  savings  by  completing  the  vectorlzation  of  the  TSTEP 
routi nes. 
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Table  VII 

PERFORMANCE  SUMMARY 


PARAMETER 

FILE1 

VAX/AP 

WALL„  CPU 

1 

WALL 

l/AX 

CPU 

PERFORMANCE 

IMPROVEMENT3 

DIRECT 

I/O  COUNT 

TIME* 

TIME 

TIME 

TIME 

WALL 

CPU 

VAX/AP 

VAX 

BM005 

57 

2.20 

98 

96.13 

1.72 

43.70 

162 

10 

BM006 

92 

19.14 

99 

97.38 

1.08 

5.09 

1862 

36 

BM008 

312 

136.12 

121 

117.64 

0.39 

0.86 

11546 

27 

NOTE:  1. 

FDAP  was  executed  with  parameter  file  BM005  and  BM006.  FDAPSEG 

was  executed  with  parameter  file  BM008. 


2.  All  times  are  in  seconds. 

3.  This  is  defined  as  VAX  time  divided  by  VAX/AP  time. 

CONFIGURATION  STATISTICS 

VAX  11/780  AP-I20B 

VMS  V2.01  APFORTRAN  V2.02 

Default  Working  Set  =  150  64K  Words  167ns.  Main  Data  Memory 

1.25  MBYTE  Memory  4K  Words  Program  Source  Memory 

No  Other  Users 
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Table  VIII 


CPU 

WALL  CLOCK 

PRIME 

SHIFT 

SECONDS 

SECONDS 

COST 

VAX  11/780 

857.3 

995.7 

$28.30 

(2  megabytes) 

CYBER -205 

23.2 

88.0 

$27.81 

(Scalar) 

CYBER -205 

17.6 

27.0 

$21.87 

(vector) 

I 


I  X 

i 


NON-PRIME 

SHIFT 

COST 

$14.15 

$14.42 

$11.34 
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